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ABSTRACT 

We report the detection of a unique CO2 ice band toward the deeply embedded, low-mass proto- 
star HOPS-68. Our spectrum, obtained with the Infrared Spectrograph onboard the Spitzer Space 
Telescope, reveals a 15.2 /xm CO2 ice bending mode profile that cannot modeled with the same ice 
structure typically found toward other protostars. We develop a modified CO2 ice profile decomposi- 
tion, including the addition of new high-quality laboratory spectra of pure, crystalline CO2 ice. Using 
this model, we find that 87-92% of the CO2 is sequestered as spherical, C02-rich mantles, while typical 
interstellar ices show evidence of irregularly-shaped, hydrogen-rich mantles. We propose that (1) the 
nearly complete absence of unprocessed ices along the line-of-sight is due to the flattened envelope 
structure of HOPS-68, which lacks cold absorbing material in its outer envelope, and possesses an 
extreme concentration of material within its inner (10 AU) envelope region and (2) an energetic event 
led to the evaporation of inner envelope ices, followed by cooling and re-condensation, explaining the 
sequestration of spherical, CO2 ice mantles in a hydrogen-poor mixture. The mechanism responsible 
for the sublimation could be either a transient accretion event or shocks in the interaction region 
between the protostellar outflow and envelope. The proposed scenario is consistent with the rarity of 
the observed CO2 ice profile, the formation of nearly pure CO2 ice, and the production of spherical 
ice mantles. HOPS-68 may therefore provide a unique window into the protostellar feedback pro- 
cess, as outflows and heating shape the physical and chemical structure of protostellar envelopes and 
molecular clouds. 

Subject headings: astrochemistry — circumstellar matter — methods: laboratory — stars: formation 
— stars: individual (HOPS-68, FIR-2) — stars: protostars 



1. INTRODUCTION 

Stars form in cores deeply embedded within dense 
molecular clouds by gravitational collapse. As an insep- 
arable part of the star formation process, released gravi- 
tational energy from the central protostar heats the sur- 
rounding infalling envelope, and angular momentum is 
shedded via powerful bipolar outflows. These feedback 
mechanisms serve to disperse protostellar envelopes, halt 
further accretion, and sculpt the physical and chemical 
properties of star-forming regions, ultimately contribut- 
ing to the evolution and destruction o f their environment 
and t he generation of turbulence (e.g. JMcKee fc Ostrikerl 
12007ft . Understanding the feedback loop between pro- 
tostars and molecular clouds is challenging; in part, 
because protostellar evolution is a stochastic process 
where violent accretion eve nts may punctuate periods 
of more moderate accretion (lHartmann fc Kenyonlll996l : 
iDunham fc VorobvovH"2012t iFischer et al.ll2012t) . Hence. 
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most observations of any given protostar only provide 
a single evolutionary snapshot, potentially leading to 
a biased interpretation. As a remedy, observers have 
long sought tracers of the chemical and physical histo- 
ries of protostars and their environments (|Charnlevlll997l : 
lArce et al"1l2007t ). One such tracer is CO2 ice. 

Since its first detection by the Infrared Astro- 
nomi cal Satellite (jD'Hendecourt fc Jourdain de Muizonl 
Il989t) . solid CO2 has proven to be an abun- 
dant and ubiquitous constituent of the interstellar 
medium. Observations of the 4.27 and 15.2 /im vi- 
brational modes toward the Galactic molecular envi- 



ronm e nts of b oth qu i escent dark cl o uds (IWhittet et al 
1998L [20071 [20091: IBergin et al.l 120051: IKnez"et~a! 
20051) and circu mstellar envelopes of low- and high- 



mass p rotostars dGerakines et al.lll999t iNummelin et all 



20011 'iBoogert et all 
Zasowski et all 120091 



20041: iPontoppidan et al.l 120081: 
ICook et all I20ll " have revealed 
a solid CO2 abundance of ^15—40% with respect to 
H2O. Chemical models indicate that gas-phase forma- 
tion of CO2 is highly inefficient at low temperature, and 
solid C0 2 is therefore widely assumed to be produced 
through chemical reactions on t he surfaces of icy grain 
mantles (jTielens fc Hagenl[l98 2i). Despite many exper- 
imental efforts to investigate the formation routes to 
solid CO2, the precise chemical pathway is still debated 
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The spectral profile of the CO2 bending mode is 
a sensitive diagnostic of the local molecular environ- 
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ment. In observations toward cold, quiescent molecu- 
lar clouds, the profile is invariably characterized by a 
relatively broad, si ngle-peaked, asym metric abso rption 
profile at 15.2 urn j Knez et all 12003 iBergin et all 120051: 
I Whittet et "all 120071 ). The quiescent profile can gener- 
ally be modeled as a linear combination of hydrogen-rich 
(H 2 0:C0 2 ) and CO-rich (CO:C0 2 ) ice mixtures, typi- 
cally dominated by the hydrogen-rich component, at low 
temperatures (T < 20 K). 

At the other extreme, along lines-of-sight toward 
highly luminous (L w 1O 4 L0), massive protostars, 
the bending mode often exhibits a do uble-peaked sub- 
struc ture near 15.15 and 15.27 fim (|Gerakines et al.l 
Il999f h The double-peaked substructure is attributed to 
Davydov splitting, a phenomenon observed in molecu- 
lar crystals h aving more than one equivalent molecule 
per unit cell (jDavvdovl 1962T). and is cha r acter is tic of 
pure, crystalline CQ 2 (lEhrenfreund et al.l 119971 119991 : 
van Broekhuizen et al.ll2006t) . The presence of pure C0 2 
is generally interpreted as the result of thermal pro- 
cessing of icy grains se en along quiescent lines-of-sight. 
IGerakines et al.l ([1999) modeled the formation of pure 
C0 2 as a segregation process of the hydrogen-rich com- 
ponent, requiring strong heating to high temperatures (T 
« 100 K, or close to the desorption temperature of the 
H 2 ice). iPontoppidan et afl (|2008| ) suggested, by mod- 
cling the presence of C0 2 ice toward low-luminosity (L 
sw 1 Lq) protostars, that pure C0 2 is also produced fol- 
lowing the thermal desorption of CO from the CO:C0 2 
mixture. This distillation process has the advantage that 
it can occur at much lower temperatures ( T = 20—30 
K). In either case, the band profile of solid C0 2 is an 
irreversible tracer of thermal processing. That is, by 
observing C0 2 ice, it can be determined whether the 
molecular cloud material has been heated significantly 
above ~20— 30 K, even if it has since re-cooled. As such, 
past accretion outbursts have recently been invoked to 
explain the formation of pure CQ 2 ice i n low-luminosity 
( L < 1 L P) ) protostars (iKim et al.ll2012l) . 

IPontoppidan et afl (|2008| ) found that all observed C0 2 
bending mode profiles toward protostellar environments 
can be decomposed, phcnomcnologically, into five unique 
components. Through the comparison with laboratory 
ice analogs, these components are ascribed to pure C0 2 
ice, C0 2 mixed with CO or H 2 ice, dilute C0 2 mixed 
with pure CO ice, and C0 2 mixed with annealed CH3OH 
and H 2 ices. Each component has a fixed profile, 
but their relative column densities vary from source- 
to-source, likely reflecting variations in the pristine-to- 
processed ice fraction along the line-of-sight. Moreover, 
this "unique component" structure appears to be a com- 
mon property of other ice bands for med near protostars 
(e.g., CO and the 6—8 11m complex; IPontoppidan et afl 
1200a IBoogert et al.ll200l . 

There exists however, at least one exception to the 
unique component model; in this paper, we demonstrate 
that the 15.2 /im C0 2 ice band profile observed toward 
the peculiar protostar, HOPS-68, cannot be modeled 
with the same C0 2 ice structure that characterizes other 
protostars. Origin ally identified as the compact millime- 
ter source FIR-2 dMezger et all 119901 ). the object is lo- 
cated in the Orion Molecular Cloud 2 region (OMC-2; 
iGatlev et al.lll974D a t an adopted di s tance of 414 ± 7 pc 
(jMenten et all 120071 ). IPoteet et all (|201lD showed that 



HOPS-68 is a moderately luminous (1.3 Lq) protostar 
with a flattened, relatively dense, infalling envelope that 
is viewed at an intermediate linc-of-sight (6 < 41°). From 
the observed spectral energy distribution (SED), they es- 
timate a mass infall rate of M = 7.6 x 10~ 6 Mq yr _1 . 

HOPS-68 is the first protostar to exhibit unambiguous 
evide nce for the presence of crystalline silicates in absorp- 
tion (jPoteet et al.l 120 ill ), indicating that the silicates in 
its protostellar envelope have undergone strong thermal 
processing (T > 1000 K). While the mechanisms respon- 
sible for such proce ssing are still not fully understood, 
IPoteet et al.l (|201lD proposed that amorphous silicates 
were annealed or vaporized within the warm inner re- 
gion of the disk and/or envelope and subsequently trans- 
ported outward by entrainment in protostellar outflows. 
Alternatively, an in situ formation by outflow-induced 
shocks may be responsible for the production of such 
material. As such, HOPS-68 might offer us a rare win- 
dow into the effects of outflows on the solid component 
of protostellar envelopes. 

In Section [21 an overview of the spectroscopic observa- 
tions, data reduction, and optical depth spectra are pre- 
sented. In Section |3l wc demonstrate that the 15.2 /im 
bending mode profile toward HOPS-68 is anomalous, and 
derive a new model for the C0 2 ice structure. Finally, in 
Section [4] we discuss possible physical scenarios that may 
explain the anomalous C0 2 ice structure toward HOPS- 
68 in the context for more typical low-mass protostars. 
Appendix [A] presents new high-quality laboratory spec- 
tra of pure, crystalline C0 2 ice that are deployed in our 
analysis. 

2. THE SPITZER-1RS SPECTRUM 
2.1. Observations and Data Reduction 

A high-resolution 10—37 /im spectrum of HOPS-68 
(a = 05 h 35 m 24?30, 5 = -05°08'30'/6 [J2000]) was ob- 
tained on 2008 Novemb e r 14 using the Spitzer Space 
Telesc ope (IWerner et al. 120041 ) Infrared Spectrograph 
(IRS; iHouck et al.1 120041 ). The observations (Spitzer 
AOR 26614016) were performed with the high-resolution 
(A/AA -600) IRS modules, Short-High (SH; 9.9-19.6 
/im) and Long-High (LH; 18.7—37.2 /im), at each of the 
two nominal nod positions, one-third of the way from the 
slit ends. The total exposure times were 120 and 240 s 
in the SH and LH modules, respectively, and were split 
between on- and off-source sky pointings. 

The spectrum was extracted from the Spitzer Science 
Center (SSC) S18.7 pipeline b asic calibrated data using 
the SMART software package (jHigdon et al.ll2004D . Per- 
manently bad and 'rogue' pixels were identified and cor- 
rected for by interpolation over nearby good pixels in 
the dispersion direction of the two-dimensional image. 
In order to remove sky emission from the high-resolution 
data we averaged the sky observations and subtracted 
their mean from the on-source data for HOPS-68. The 
sky subtracted two-dimensional images were averaged by 
nod position, and extracted using a full slit extraction. 
Spectra of Markarian 231 (LH) and a Lac (SH) were pro- 
duced in the same manner, and relative spectral response 
functions (RSRFs) were produced by dividing a template 
by these spectra. Each nod of the HOPS-68 spectrum 
was then multiplied by these RSRFs to calibrate the flux 
density scale and the nods were averaged to obtain a final 
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Figure 1. Top panel: Spitzer-TRS 13—19 fim SH spectrum of 
HOPS-68, shown with statistical uncertainties (gray lines) and the 
adopted continuum (thick gray line). The broad ab sorption near 
16.1 lira has been ascribed to forsterite (Mg2Si04; IPoteet et al.l 
I201U ). and the narrow lines at 16.0 and 17.03 (im are background- 
subtraction artifacts; the latter resulting from H2 5(1) emission 
near the source. Bottom panel: derived optical depth spectrum 
of HOPS-68 (black line), centered near the strong 15.2 /an CO2 
ice V2 bending mode. The sharp absorption feature at 14.97 (im 
corresponds to the blended Q-branch ro- vibrational lines of gaseous 
CO2. Insert: close-up of the 14.5—16 fim spectrum of HOPS-68 
(black line), with scaled optical depth profiles from the protostars: 
S140 IRS 1 (blue line; x9.0), RNO 91 (green line; x3.4), and IRAS 
03254 (magenta line; x4.5) shown for comparison. 



spectrum. Finally, the SH and LH spectra were scaled 
to the flux density of the second-or der of the Long-Low 
spectrum (LL2; 14.2-20.4 fim) from lPoteet et al.l (|2011l ). 
The resultant SH spectrum is shown in the top panel of 
Figure [TJ 

2.2. Optical Depth Spectrum 

We convert the CO2 flux density to an optical depth 
spectrum, using: r(A) = - ln(f ° bs /F™ nt ), where F° bs 
and F%° nt are the observed and continuum flux densi- 
ties, respectively. The continuum is constructed by fit- 
ting a third-order polynomial to the observed flux density 
within the wavelength ranges of 13.0—13.3, 14.6—14.7, 
and 18.2—19.5 fim; the 13.5—14.2 fim range is avoided 
due to the rare gas-phase absorption of C2H2 (13.71 fim) 
and HCN (14.05 fim). The polynomial is then combined 
with a Gaussian in the wavenumber domain, with a cen- 
ter at v = 608 cm" 1 and a Full Width at Half Maximum 
(FWHM) of Ai> = 73 cm" 1 , to sim ulate the blue wing 
of th e 18 fim silicate bending mode ([Pontoppidan et all 
12008ft . 



3. AN ANOMALOUS C0 2 ICE BAND 

The derived CO2 ice optical depth spectrum of HOPS- 
68 is presented in the bottom panel of Figure [TJ The 
profile exhibits prominent double-peaked substructure at 
15.05 an d 15.24 fim, characteristic of pure, crystallin e 
C0 2 ice dHudgins et all 119931: lEhrenfreund et al.lll997ft . 
The C0 2 bending mode profiles toward the well-studied 
protostars: S140 IRS 1, RNO 91, and IRAS 03254 are 
shown for comparison in the insert of Figure [TJ illustrat- 
ing that the HOPS-68 pr o file is fundamentally d ifferent 
(jPontoppidan et al.l [2008t iGerakines et al.l 11999ft . The 
double-peaked substructure toward these low- and high- 
mass protostars occurs near 15.1 and 15.23 fim and is less 
pr ominent in comparison ; the dip-to-pcak ratios, defined 
by IZas owski et al. (2009|) as the local minimum to local 
maximum ratio of the blue peak, is 1.6—3.4 times less 
than that of HOPS-68. Secondly, toward these protostars 
a third peak or broad shoulder is often observed near 15.4 
fim, but is only weakly detected in the SH spectrum of 
HOPS-68. Finally, the SH spectrum of HOPS-68 reveals 
the presence of both gaseous 14.97 fim CO2 absorption 
and a blue absorption component or wing shortward of 
15.05 fim. In comparison, gas-phase C0 2 absorption is 
weakly observed in the spectrum of the massive protostar 
S140 IRS 1; however, the spectrum displays no evidence 
for a blue wing. 

3.1. Unique Component Model 

We first demonstrate that the unique component de- 
composition, which generally fits the CO2 bending mode 
profiles in both low- a nd high-mass proto s tars, fails 
for HOPS-68. Following IPontoppidan et al.l (|2008[) . we 
use the following five components, corrected for shape- 
effects using a continuous distribution of ellipsoids (CDE; 
iBohren &: Huffman] I1983D model: hydrogen-rich C0 2 
(H 2 0:C0 2 = 100:14; T = 10 K), CO-rich C0 2 (CO:C0 2 
= 100:70, 100:26 or C0 2 :CO = 112:100; T = 10 K), pure 
C0 2 (T — 15 K), and dilute C0 2 (CO:C0 2 = 100:4; T 
= 10 K). For the pure C0 2 component, we adopted the 
shape-corrected, high-resolution spectra from this work 
(see Appendix To match the resolving power of the 
SH module, the high-resolution spectra are convolved 
with a Gaussian, having a FWHM of AA « 0.03 fim, 
at each wavelength element. The fifth component, often 
referred to in the literature as the CO2 shoulder at 15.4 
fim, has been identified as an intermolecular interaction 
between C Q 2 and CH 3 OH in annea l ed HaO^HaOHrCO a 
mixtures (|Ehrenfreund et alj[l99l iDartois et all 11999ft . 
and is empirically modeled using a superposition of two 
Gaussians in the wavenumber domain. 

We em ploy a nonlinear least-squares fitting routine 
(MPFIT; lMarkwardtl2009ft to determine the best-fit scal- 
ing of each component to the observed CO2 bending 
mode of HOPS-68: 

5 

r(A) = Y J ociT i (\), (1) 

?:=i 

where Ti(X) is the normalized optical depth of an individ- 
ual laboratory ice component and on is its corresponding 
scaling parameter. The fit is performed over the wave- 
length ranges of 14.75—14.95 and 15.0—15.6 fim to avoid 
the gas-phase CO2 absorption at 14.97 fim and the blue 
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Figure 2. Spitzer-YRS SH optical depth spectrum o f HOPS-68 (black li n e) an d best-fit CO2 ice models. Left panel: the best-fit unique 
component CDE grain shape model (blue line) from Pontoppidan et al. (2008). The model includes contributions from the pure CO2 
component (T = 15 K; green line), the hydrogen-poor C02:CO = 112:100 mixture (T = 10 K; viol e t line ), and the 15.4 (im shoulder 
component (orange line). Right panel: the best-fit annealed ice model adopted from IGe rakines ct al. (1999). The model consists of the 
annealed H20:CH30H:C02 = 1:0.9:1 mixture (T = 125 K; green line) and the CDE shape-corrected laboratory spectrum of the hydrogen- 
rich H20:C02 = 100:14 mixture (T = 10 K; red line). The fitting procedures are constrained to wavelengths shortward of 15.6 fim to 
avoid the blue wing of the 16.1 fim forsterite feature. 



wing of the 16.1 /xm forsterite feature. 

A comparison between the best-fit model and the 
observed CO2 bending mode of HOPS-68 is shown 
in the left panel of Figure [21 and clearly indicates 
that the ice composition is poorly simulated (jq fa 
391) by the unique co mponent decomposition from 
iPontoppidan et al.l (|2008t ). The best- fit model results in 
a three-component ice mantle that is dominated by pure 
CO2, with contributions from the CO-rich C02:CO = 
112:100 mixture and the 15.4 /im shoulder component. 
The hydrogen-rich H20:C02 = 100:14 mixture, which 
generally dominates the bendi ng mode profiles toward 
low- and high-mass pr o tostars ([Pontoppidan et al.ll2008t 
Gerakines et all 119991: lOliveira et all 120091: i Seale et al.l 
20111: lAn et alJl20lH iKim et al.ll20T2n . is entirely absent 
from the model. The red wing of the observed profile is 
adequately matched by the pure CO2 component, how- 
ever the blue peak (15.14 /xm) of the pure CO2 spec- 
trum does not coincide with the observed peak position of 
15.05 /xm. Moreover, we find that the CO-rich C02:CO 
= 112:100 mixture is too broad and redshifted to account 
for the observed blue wing of the bending mode; decreas- 
ing the CO2 fraction of this mixture results in a narrow 
profile that becomes further redshifted. 

3.2. Annealed H 2 0:CH 3 OH:C0 2 Model 

We now consider the approach from iGerakines et al.1 
(jl999f ) , who first interpreted the CO2 bending mode pro- 
files toward high-mass protostars using a combination 
of a hydrogen-rich mixture and an annealed CH3OH- 
rich mixture. Similarly to ISeale et"aT] (|2011l) . a two- 
component model was assembled using the CDE shape- 
corrected H 2 0:C0 2 = 100:14 (T = 10 K) mixture and 
an annealed H20:CH30H:C02 mixture fr om the large 
library of laboratory spectra described in IWhite et al.l 
(|2009t) . However, since the blue peak of the observed 
bending mode is positioned near 15.05 /im, we con- 
sider only the bluest {y > 662 crn^ 1 or A < 15.1 /xm) 
H 2 0:CH 3 0H:C02 spectra for simplicity. The annealed 



mixtures are describ ed in Table [U For cons i stency with 
other studies (e.g., IGerakines et al.l 119991 : iCook et al.l 
120111) . grain shape corrections were not applied to the 
absorption spectra of annealed mixtures; the segregation 
of CH3OH and CO2 in laboratory ice mixtures is thought 
to result in a highly inhomogeneous str ucture that is no 
longer well represented by th in films ([Gerakines et all 
[i99aiEhrenfreund et al J [19991) . 

The best-fit model is presented in the right panel 
of Figure [H In contrast to past studies (e.g., 
IGerakines etaTI 119991: ISeale et al.l 1201 ID . we find that 
combining the CDE shape-corrected hydrogen-rich 
H 2 0:C0 2 = 100:14 (T = 10 K) mixture with an an- 
nealed H 2 0:CH 3 OH:C02 = 1:0.9:1 (T = 125 K) mix- 
ture yields a poor goodness of fit (xt ~ 235). In par- 
ticular, we find that the hydrogen-rich mixture, which 
is known to consistently dominate the bending mode, 
is simply not utilized by the nonlinear least-squares fit- 
ting routine due to its broad redshifted profile. Instead, 
the fitting routine attempts to constrain the ice compo- 
sition with a single annealed mixture. However, even 
for the bluest annealed mixtures, we find that the po- 
sition of the blue peak does not coincide with the ob- 
served peak position of ~15.05 (im. This result demon- 
strates that the observed CO2 bending mode cannot be 
modeled by combining a hydrogen-rich mixture with an 
annealed H20:CH30H:C02 mixture, or by a single an- 
nealed H20:CH30H:C02 mixture. Moreover, this is con- 
sistent with the fact that the unique profile decomposi- 
tion is capable of simulati ng the same lines-of-sigh t as 
the annealed ice mixtures (jPontoppidan et al.|[2008l) . 

3.3. Multi- Component Homogeneous Spheres Model 

Having found that the traditional methods for ana- 
lyzing the CO2 ice band profile fail for HOPS-68, we 
now investigate a broader set of possibilities. In par- 
ticular, we must include a component that can provide 
additional absorption in the blue wing of the profile. To 
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Figure 3. Spitzer-IRS SH optical depth spectrum of HOPS-68 (black line) and best-fit CO2 ice models (blue lines). Left panel: the best-fit 
multi-component homogeneous spheres model. Individual CO2 ice components include: hydrogen-poor CO2 C02:H20 = 6:1 mixture (T 
= 55 K; violet line), pure CO2 (T = 15 K; green line), hydrogen-rich H20:C02 = 100:14 mixture (T = 10 K; red line), and the 15.4 fim 
shoulder (orange line). Right panel: the best-fit multi-component CDE + spheres model. A line-of-sight consisting of two grain shape 
populations was simulated by combining the CDE shape-corrected laboratory spectra of pure CO2 ice (T = 75 K; green line) and the 
hydrogen-rich H20:C02 = 100:14 mixture (T = 10 K; red line), with the homogeneous C02:H20 = 6:1 spheres mixture (T = 55 K; violet 
line). 



this end, we consider additional CQ 2 -bea r ing, h ydrogen- 
poor ice analogs from lEhrenfreund et al.l (|1997l) and ex- 
plore the effect of grain shape by using homogeneous 
spheres in the Rayleigh limit. The hydrogen-poor mix- 
tures, which are summarized in Table [1] consist of two or 
more molecular constituents and span a laboratory tem- 
perature range of T = 10—80 K. Grain shape-corrected 
spectra for a variety of particle s hape distributions are 
presented by lEhrenfreund et all (|1997f ). For relatively 
dilute C0 2 (e.g., H 2 0:C0 2 = 100:14 and CO:C0 2 = 
100:26), grain shape effects are generally weak, and only 
minor differences in the bending mode profile exist be- 
tween the CDE and spherical grain shape corrections. 
Conversely, for C0 2 -rich mixtures (e.g., pure C0 2 and 
C0 2 :H 2 = 6:1), spherical grain shape corrections often 
result in a narrow blueshiftcd profile that is qualitatively 
consistent with the observed blue wing. 

Following a similar methodology as before, a spec- 
tral decomposition of the observed bending mode 
was performed usi n g the multi-component model from 
iPontoppidan et all (|2008[ ). However, for this model, 
the choice of the hydrogen-poor C0 2 component was 
not restricted to t h e thr ee CO-rich mixtures used by 
IPontoppidan et all (|2008[ ). but instead selected from 
the aforementioned suite of C0 2 -bearing ice analogs by 
lEhrenfreund et all (|1997| ). The laboratory temperature 
of the pure C0 2 component was also permitted to vary 
from T = 15—75 K. Furthermore, we assume that the 
line-of-sight consists exclusively of small homogeneous 
spheres. 

The best-fit homogeneous spheres model is presented in 
the left panel of Figure [3] A four-component model, con- 
sisting largely of pure C0 2 ( T = 15 K) and the hydrogen- 
poor C0 2 :H 2 = 6:1 (T = 55 K) mixture, was found to 
provide the best fit (xl ~ 19) to the observed bend- 
ing mode. As qualitatively expected, we find that the 
observed blue wing is adequately simulated by spheri- 
cal grain shape-corrected laboratory spectra of C0 2 -rich 



mixtures. Furthermore, the observed red wing is well 
matched by a combination of pure C0 2 , the hydrogen- 
rich H 2 0:C0 2 = 100:14 (T = 10 K) mixture, and the 
15.4 yLtm shoulder component. However, due to the rela- 
tively low amplitude of the pure C0 2 red peak, the opti- 
cal depth of the 15.24 /mi peak is slightly under-predicted 
by the best-fit model. Nonetheless, we conclude that an 
exclusive homogenous spheres model, consisting predom- 
inately of the hydrogen-poor C0 2 :H 2 = 6:1 (T = 55 
K) mixture, yields an overall excellent fit to the observed 
C0 2 bending mode profile. 

3.4. Homogeneous Spheres + CDE Model 

The grain shapes of CO and C0 2 ice mantles in 
the circumstellar envelopes of protostars have been 
constrai ned using multiple ice bands of the same 
species dBooeert et all 120021 : IPontoppidan et al.l 120031 : 
iGerakines et al .111 9991 ). A simultaneous decomposition of 
these bands indicates that their observed profiles are con- 
sistent with those of small irregularly-shaped ice mantles, 
simulated by CDE grain shape models. For this reason, 
invoking a model consisting exclusively of homogeneous 
spheres may seem unjustified. However, wc have demon- 
strated that utilizing a CDE grain shape model to sim- 
ulate the C0 2 bending mode of HOPS-68 is largely un- 
satisfactory. Even when the choice of the hydrogen-poor 
component is no t constrained to the CO- rich CO:C0 2 
mixtures used bv lPontoppidan et al.l (|2008f ). we find that 
the goodness of fit is generally poor (xt ~ 90—680). Al- 
though our analysis is restricted to laboratory experi- 
ments for which grain shape corrections have been ap- 
plied, this result implies that the observed blue wing 
is inconsistent with all available CDE shape-corrected 
hydrogen-poor ice analogs. In contrast, we have shown 
that small homogeneous spheres of the hydrogen-poor 
C0 2 :H 2 = 6:1 (T = 55 K) ice mixture can adequately 
fit the blue wing of the observed bending mode. 

To further examine the robustness of this result, wc 
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Table 1 

Summary of Laboratory Absorption Spectra 



Ice 

Composition 



Mixture 
Ratio 



T Resolution 
(K) (cm" 1 ) 



Hydrogen-rich Ice (Ehrcnfrcund et al. 1997) 



H 2 0:C0 2 



100:14 



10 



Hydrogen-poor Ices fEhrcnfrcundetal. 1997) 

C():C()., 

CO:C0 2 

CO:C0 2 

CO:C0 2 

CO:C0 2 

CO:C0 2 

CO:C0 2 

CO:C0 2 

CO:C0 2 

CO:C0 2 

CO:C0 2 

CO:C0 2 

CO:C0 2 

CO:N 2 :C0 2 

CO:N 2 :C0 2 

CO:0 2 :C0 2 

CO:0 2 :C0 2 

CO:0 2 :C0 2 

CO:0 2 :C0 2 

CO:0 2 :C0 2 

CO:0 2 :C0 2 

CO:0 2 :C0 2 

CO:0 2 :C0 2 

CO:0 2 :C0 2 

CO:0 2 :C0 2 

CO:0 2 :C0 2 

CO:0 2 :C0 2 

CO:0 2 :C0 2 

CO:0 2 :C0 2 

CO:0 2 :C0 2 

CO:0 2 :C0 2 

CO:0 2 :N 2 :C0 2 

CO:0 2 :N 2 :C0 2 

CO:0 2 :N 2 :C0 2 :H 2 . 
CO:0 2 :N 2 :C0 2 :H 2 . 
CO:0 2 :N 2 :C0 2 :H 2 . 
CQ 2 :CO 



C0 2 :CO . 
C0 2 :H 2 
C0 2 :H 2 
C0 2 :H 2 
C0 2 :H 2 
C0 2 :H 2 
C0 2 :H 2 
C0 2 :H 2 
C0 2 :H 2 
C0 2 :H 2 
C0 2 :H 2 
C0 2 :0 2 . 



100:4 


10 


1 


100:4 


30 


1 


100:8 


10 


1 


100:8 


30 


1 


100:16 


10 


1 


100:16 


30 


1 


100:21 


10 


1 


100:21 


30 


1 


100:23 


10 


1 


100:23 


30 


1 


100:26 


10 


1 


100:26 


30 


1 


100:70 


10 


1 


100:50:20 


10 


1 


100:50:20 


30 


1 


100:10:23 


10 


1 


100:10:23 


30 


1 


100:11:20 


10 


1 


100:11:20 


30 


1 


100:20:11 


10 


1 


100:20:11 


30 


1 


100:50:4 


10 


1 


100:50:4 


30 


1 


100:50:8 


10 


1 


100:50:16 


10 


1 


100:50:16 


30 


1 


100:50:21 


10 


1 


100:50:21 


30 


1 


100:50:32 


10 


1 


100:54:10 


10 


1 


100:54:10 


30 


1 


100:50:25:32 


10 


1 


100:50:25:32 


30 


1 


25:25:10:13:1 


10 


1 


50:35:15:3:1 


10 


1 


50:35:15:3:1 


30 


1 


112:100:1 


10 


1 


112:100:1 


45 


1 


6:1 


10 


1 


6:1 


12 


1 


6:1 


15 


1 


6:1 


50 


1 


6:1 


55 


1 


6:1 


75 


1 


10:1 


10 


1 


10:1 


80 


1 


100:1 


10 


1 


100:1 


30 


1 


1:1 


10 


1 



Pure C0 2 Ices (This Work) 



Pure C0 2 




15 


0.1 


Pure 0O 2 




30 


0.1 


Pure C0 2 




45 


0.1 


Pure 0O 2 




60 


0.1 


Pure CQ 2 




75 


0.1 


Annealed Ices (White et al. 2009) 


H 2 0:CH 3 OH:C0 2 


1:0.9:1 


125 


1 


H 2 0:CH 3 OH:C0 2 


1:0.9:1 


130 


1 


H 2 0:CH 3 OH:C0 2 


1:1.7:1 


110 


1 


H 2 0:CH 3 OH:C0 2 


1:1.7:1 


115 


1 


H 2 0:CH 3 OH:C0 2 


1:1.7:1 


120 


1 



next consider a modeled linc-of-sight that is populated by 
both spherical and irregularly-shaped grains. Following a 
similar approach as in Section [3.31 a CO2 ice model was 
constructed by combining the previously utilized CDE 
shape-corrected laboratory spectra (hydrogen-rich CO2, 
pure CO2, and dilute CO2), a spherical grain shape- 
corrected hydrogen-poor laboratory spectrum, and the 
empirical 15.4 /im CO2 shoulder component. The best- 
fit model is shown in the right panel of Figure [3] Indeed, 
we find that a four-component model, composed largely 
of the same homogeneous CC^H^O = 6:1 (T = 55 K) 
spheres mixture from Section [3.31 is required for a satis- 
factory fit (xl « 24). The red wing of the observed bend- 
ing mode is well matched by the combination of the CDE 
shape-corrected spectra of pure CO2 ( T — 75 K) and the 
hydrogen-rich H 2 0:C0 2 = 100:14 (T = 10 K) mixture, 
with the 15.4 /im shoulder component. In comparison 
to the model results from Section l3~3[ we conclude that 
a modeled line-of-sight composed of two separate grain 
shape populations (spherical and irregularly-shaped par- 
ticles) provides an equally good fit. This result demon- 
strates that the CO2 bending mode profile toward HOPS- 
68 is largely governed by the grain shape and chemical 
composition of the hydrogen-poor CO2 component. 

3.5. CO2 Ice Column Density 

The column density along the line-of-sight for each 
best-fit CO2 component was estimated from Equation 
IA2[ adopting the pure CO2 bending mo de band strength 
of A = 1.1 x 10~ 17 cm molecule -1 from lGerakines et al.l 
(|1995h . The optical depth spectra were integrated 
over the wavenumber range of 680.3—613.5 cm -1 (A = 
14.7—16.3 /im) to include the long wavelength wing of 
the observed profile. 

To estimate the systematic contribution from the base- 
line subtraction, we first computed the standard error, 
a, between the best-fit continuum and the SH spectrum. 
The best-fit continuum was then offset by ±1 a to pro- 
duce new optical depth spectra and column density esti- 
mates. To quantify the errors due to the baseline uncer- 
tainty, we then calculated the mean absolute difference 
between the best-fit column densities and the column 
densities obtained by offsetting the best-fit continuum. 
These errors were then added in quadrature with the 
best-fit errors, which are returned by the least-squares 
minimization technique, to generate the final CO2 ice 
column density uncertainties. 

The column densities and relative abundances for the 
best-fit CO2 ice models are tabulated in Table [2] The 
modified spectral decomposition analysis gives us the fol- 
lowing three observablcs: 1) Among the available labora- 
tory ice analogs, homogeneous spheres of the hydrogen- 
poor C02:H20 = 6:1 (T = 55 K) mixture are required 
to match the blue wing of the observed CO2 bending 
mode. This component constitutes 65—74% of the total 
CO2 ice column density along the line-of-sight. 2) The 
red wing of the observed profile can be modeled by com- 
bining pure CO2 ice and the 15.4 /im CO2 shoulder com- 
ponent. The pure CO2 component constitutes 18—22% 
of the total CO2 ice column density, while only a minor 
(3—6%) contribution is attributed to the 15.4 /im CO2 
shoulder component. 3) Only 2—10% of the total CO2 
ice column density is composed of the hydrogen-rich com- 
ponent, while no evidence for the CO-rich component is 
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Table 2 

Column Densities and Abundances for Best-fit CO2 Ice Models 



Grain Shape 
Model 


AT(C0 2 ) 
(10 17 cm" 2 ) 


H-rich CO-rich Pure Shoulder 

(%) (%) (%) (%) 


Dilute H-poor 

(%) (%) 


(K) 


xi 


Multi-component Models 


CDE 

Spheres 

CDE+Spheres 


19.83 ± 0.49 
24.66 ± 0.59 
. 24.53 ± 0.59 


27 ± 3 65 ± 3 8 ± 12 
10 ± 13 ■ ■ ■ 22 ± 3 3 ± 11 
2 ± 52 ■■■ 18 ±3 6 ± 8 


65 ± 3 
74 ± 3 


15 
15 

75 


391.1 
18.6 
24.3 


Median Abundance Toward Low-mass Protostars (Pontoopidan et al. 2008") a 


CDE 




roll 20H 6 1 ? 31 


<0.ll 


15 





Note. — Uncertainties are based on statistical errors in the Spitzer-S}! spectrum and systematic errors 
in the continuum determination. Grain shape models from this work include: irregular grains (CDE), pure 
homogeneous spheres (spheres), and two separate grain shape populations (CDE+Spheres). For the CDE 
and spheres models, all CO2 components consist of irregularly-shaped particles and homogeneous spheres, 
respectively. For the CDE+sphcres model, all components are representative of irregularly-shaped particles, 
except for the hydrogen-poor CO2 component which is representative of homogeneous spheres. 
a Median abundances include upper and lower quartile values. 



found within the observed CO2 bending mode. These 
fundamental differences of HOPS-68, as compared to the 
ice structure toward typical low-mass protostars, are il- 
lustrated in Table [2] 

To derive the CO2 ice abundance, the H 2 ice column 
density toward HOPS-68 was estimated from the 12 fim 
libration mode of the amo rphous silicate-subtr acted op- 
tical depth spectrum from iPoteet et al.l (|2011l) . Adopt- 
ing an intrinsic band strength of A = 2.9 x 1CP 17 
cm molecule" 1 for c rystalline H 2 ice (T = 140 K; 
iMastrapa et al.ll2009f ). we calculate a H 2 ice column 
density of iV(H 2 0) = 88.2 ± 0.7 x 10 17 cm" 2 . For 
the best-fit models from Sections 13.31 and 13. 4[ we cal- 
culate an average C0 2 to H 2 ice column density ratio 
of #(C0 2 )/iV(H 2 0) = 0.28 ± 0.03. This is consistent 
with abundances measured toward l ow-mass protostars 
HV(C0 2 )/A r (H 2 0) = 0.32 ± 0.02; IPontoppidaneTall 
mim . but is significantly greater than those observed for 
high-mass protostars ( N (CO?) / N (K?.Q) = 0.17 ± 0.03; 
iGerakines et al.lll999t) . 

4. DISCUSSION AND CONCLUSIONS 

4.1. The Rarity of Highly Processed Ice 

The Spitzer-IRS SH spectrum toward the low-mass 
protostar HOPS-68 reveals a rare 15.2 fim C0 2 bend- 
ing mode profile, exhibiting a prominent double-peaked 
substructure (near 15.05 and 15.24 fim) and a strong 
blue absorption wing. Our analysis demonstrates that 
the observed bending mode cannot be modeled using 
the un ique component models from iPontoppidan et al.l 
(|2008[ ) or by combining a hydrogen-ri ch mixture with 
an an nealed H 2 0:CH3QH:CQ 2 mixtu re (Gerakinc s et al.l 
119991) . However, when the IPontoppidan et al.l (|200"8f ) 
models are modified to include spherical grains and a 
broader selection of hydrogen-poor ice mixtures, an over- 
all excellent fit to the observed C0 2 bending mode is 
obtained. 

Our results indicate that most (87—92%) of the C0 2 
along the line-of-sight is present in ices where C0 2 is 
the dominant constituent. In contrast, the majority 
of the C0 2 ice observed toward most low- and high- 
mass protostars is found in environments dominated by 
H 2 and CO. For example, toward low-mass protostars, 



IPontoppidan et all (|2008l) showed that generally ^70% 
of the total C0 2 ice column density is composed of the 
hydrogen-rich component, while ~20% is attributed to 
the CO-rich component. Furthermore, the median abun- 
dance of pure C0 2 ice to the total C0 2 ice column den- 
sity is found to be ~6%. 

The nearly complete lack of unprocessed ice mantles 
along the line-of-sight toward HOPS-68 is unexpected. 
For example, toward the prominent low-mass sources 
IRAS 03254 and RNO 91, comp arable fractions of 
pure C Q 2 (17—23%) are reported by IPontoppidan et al.l 
i|2008D . However, a substantial fraction (70-75%) 
of their total C0 2 ice column density is attributed 
to unprocessed ice mixtures (CO:C0 2 and H 2 0:C0 2 ), 
which lie outside the distillation and segregation radii 
(jOberg et al.l l20iTl ). In this context, the ice mantles 
toward HOPS-68 are uniquely different from those ob- 
served toward other protostars, indicating that the ob- 
served icy grains have undergone a higher degree of ther- 
mal processing. 

4.2. The Origin of Processed Ice Toward HOPS-68 

Nearly 100 C0 2 ice spectra have been described in the 
literature, and the observed profile toward HOPS-68 is 
unique among them. It is therefore unlikely to find more 
than a few percent of HOPS-68 ice analogs in a larger 
sample of protostars. Is this rarity due to HOPS-68 
being fundamentally different than other protostars, or 
does the observed line-of-sight reveal processed ices that 
are present in other protostars, but which are typically 
obscured by their cold outer envelope? 

4.2.1. The Lack of Unprocessed Ice 

A peculiar spectral property of HOPS-68 is the combi- 
nation of a relatively flat 8 —70 fim SEP and strong sili- 
cate absorption at 9.7 fim (jPoteet et alj|201ll Figure 3) 
To rep roduce these observed characteristics . IPoteet et all 
(| 2 1 If) employed a radiative transfer model based on the 
collapse of a n isothermal sheet ini t ially in hydrostatic 
equilibrium (|Hartmann et al.l 11994 |1996[) . They con- 
clude that HOPS-68 possesses a highly flattened (r] — 
2.0) protostellar envelope, as described by the degree 
of asphericity (rj = R max /H, where iJ max is the outer 
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envelope radius and H is the scale height o f the initial 
sheet). In comparison, iFurlan et al.l (120081 ) found that 
most Taurus protostars do not require an initially flat- 
tened density distribution, but instead can be modeled 
by th e collapse of an initia lly spherically symmetric cloud 
core (jTerebev et al.l [l984[) . Unlike spherical envelopes, 
an implication of the highly flattened structure is that 
cold material in the outer envelope is not expected to 
be present along the observed line-of-sight. Indeed, we 
estimate that <15% of the gas has temperatures below 
T = 30 K. Thus, the absence of the volatile CO-rich 
CO2 component and the low abundance of the hydrogen- 
poor CO2 component may be explained by the lack of 
cold, unprocessed dust along the line-of-sight. This re- 
sult suggests that highly processed ices may be present 
in the inner regions of other protostcllar envelopes, but 
arc likely obscured or diluted by the large, cold envelopes 
of unprocessed material typically surrounding them. 

Furthermore, to si multaneously s imula te the strong sil- 
icate absorption, the lPoteet et al.l (|2011[ ) model requires 
infalling material to be deposited within 0.5 AU of the 
central protostar, as set by the centrifugal radius (i? c ). 
In contrast, R c values from 10 to 300 AU, with a median 
of 60 AU, were found for Taurus protostars (IFurlan et al.l 
200$. For a centrifugal radius of R c = 0.5 AU, we esti- 
mate that more than 50% of the infalling material along 
the line-of-sight is concentrated within 10 AU of the cen- 
tral protostar, where the potential for thermal processing 
of icy grains by radiation and outflow-induced shocks is 
the strongest. 

Resolved images of protostellar envelopes show a wide 
range of envelope morphologies, from highly flattened to 
relative sph eroidal structures, as we l l as a great deal of 
complexity (jTobin et alJl2010l 120 111 120121 ). The diver- 
sity of envelope morphologies may be due to formation 
of protostars in turbulent and often filamentary molecu - 
lar clouds (e.g., lAndre et al.ll2010t iMolinari et al.ll2010D . 
Although the envelope of HOPS-68 has not yet been re- 
solved, SED modeling requires a highly flattened enve- 
lope which may be uncommon within the diversity of en- 
velope morphologies typically found in molecular clouds. 
This flattened envelope, when observed from an inter- 
mediate inclination, results in a line-of-sight that is less 
obscured by unprocessed material in the cold, outer en- 
velope. 

4.2.2. The Formation of Spherical C02-rich Mantles 

Although a flattened envelope structure may explain 
the absence of cold primordial ices seen toward other 
protostars, it does not directly explain why a large frac- 
tion of the CO2 ice seen toward HOPS-68 is sequestered 
as spherical, C02-rich mantles, while typical interstel- 
lar ices show evidence of irregularly-shaped, hydrogen- 
rich mantles. That is, in the absence of unprocessed 
ices along the line-of-sight, the observed CO2 ice profile 
should be consistent with irregularly-shaped C02-rich ice 
mantles. Irregularly-shaped ice mantles are generally un- 
derstood to be the product of a slow aggregation of icy 
grain monomers in dense molecular clouds prior to pr oto- 
stellar collapse (jOssenkopj 119931 : lOrmel et al.M2009Q . In 
the case of HOPS-68, we propose that the spherical ice 
mantles may have formed as the volatiles rapidly freeze 
out in dense gas, following an energetic, but transient 
event that sublimated any primordial icy grain popula- 



tion within its inner envelope region. 

Radiative processing by episodic accretion events have 
recently been proposed by iKim et al.l (|2012l ) to be a vi- 
able mechanism for producing pure CO2 ice, through 
the distillation process, toward low-luminosity protostars 
(L < 1 L p)). Using the uniq u e component decomposi - 
tion from iPontoppidan et al] (|2008l) . IKim et al.l (|2012f ) 
report a median pure CO2 ice abundance of 15% to- 
ward six protostars; this value is six times less than the 
processed ice abundance found toward HOPS-68. The 
detection of pure CO2 ice toward low-luminosity proto- 
stars indicates that a transient phase of high er luminos- 
ity m ust have existed some time in their past (IKim et al] 
12011) . For the case of HOPS-68, we propose that heat- 
ing beyon d the sublimation temperature of H2O ice (T 
« 110 K; iFraser et al.l l200ll ) will evaporate ice mantles 
close to the central protostar during the high-luminosity 
phase, resulting in the separation of icy grain conglom- 
erates. As the protostar later returns to its quiescent, 
low-luminosity phase, gas temperatures in the envelope 
decrease and H 2 and CO2 rapidly {t < 10 3 yr) condense 
onto grain s urfaces in the or der of their condensation 
temperature (|OssenkopJ 119931 ) . provided that the gas is 
sufficiently dense (riH 2 > 10 4 cm -3 ). This mechanism 
naturally allows for the formation of nearly pure CO2 ice 
mantles sequestered on spherical grains. However, this 
scenario requires that the accretion outburst must have 
occurred recently to avoid the reproduction of irregular 
grains by renewed grain aggregation. 

Evidence for transient heating is often presumed from 
the variable nature of the accretion luminosity. Photo- 
metric monitoring of HOPS-68, with the Herschel Space 
Observatory, reveal a ^21% variation in the 70 /im flux 
over a two week period; this variability is most likely 
driven by variatio ns in the mass acc retion rate onto the 
central protostar (TBillot et al.l 120121) . While lar ger in- 
creases in luminosity, over longer timescales, are certainly 
needed to heat dust grains to the sublimation temper- 
ature of H2O ice, the observed variability indicates a 
rapidly changing accretion rate for HOPS-68, and hence, 
supports the possibility of a recent (t <C 10 6 yr) accretion 
outburst that raised the luminosity significantly above 
the current observed value. 

In order to quantify the luminosity rise needed to pro- 
duced the proposed sublimation o f icy grains, we in - 
crease the total luminosity of the iPoteet et al.l (|2011[ ) 
SED model by a factor of 10, 50, and 100. Adopting 
a CO2 ice desorption temperatu re of T = 45 K, from 
an amorphous H2O ice surface (|Noble et al.l 120121 ). we 
assume that CO2 ice is present in the envelope of HOPS- 
68 at distances greater than r « 90 AU. These incre- 
mental steps in luminosity result in increased dust tem- 
peratures of T = 70, 97, and 113 K at r « 90 AU, re- 
spectively. Thus, for HOPS-68, we propose that a fac- 
tor of 100 increase in luminosity (^100 L@) is required 
to thermally desorb CO2 and H2O ices from the sur- 
face of dust grains. Although the required luminosity 
is higher than the typical luminosit ies of protostars in 
Orion (~1 L© : iKrvukova et ~aT]|2012f) . luminosities in ex- 
cess of 100 Lq are observed toward yo ung, low-mass stars 
unde r going FU Orions- l ike ou tbursts (jReipurth fc Aspiiil 
120101 ) . IKrvukova et al.l (|2012f ) finds five Orion protostars 
with luminosities exceeding 100 L©, including two can- 
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Figure 4. Model predictions for the (£3) 12 CC>2 and (S3) 13 C02 
stretching modes at 4.27 and 4.38 (im, respectively. The predicted 
CO2 profiles (blue lines) are constructed from th e be st-fit multi- 
component homogeneous spheres model (Section 13.31 top panel) 
and the mu lti-component CDE + homogeneous spheres model 
(Section 13.41 bottom panel). The individual ice components in- 
clude: hydrogen-rich CO2 (H 2 0:C0 2 = 100:14 at T = 10 K; red 
lines), hydrogen-poor CO2 (C0 2 :H 2 = 6:1 at T = 55 K; violet 
lines), and pure CO2 (T = 15 and 75 K; green lines). 

didatc FU Orionis-likc objects. Furthermore, models de- 
veloped to explain the typical low luminosities of pro- 
tostars require that most low-mass protostars undergo 
FU Orionis-like outbursts in order to accrete the re- 
quired amount of mass dur ing the protostellar lifetime 
(Dunham fe VorobvovllMl ) . 

Alternatively, mechanical proce ssing by shocks, asso - 
ciated with the observed outflow ([Williams et al.|[2003l ). 
will also sublimate ice mantles in the outflow working 
surface along the cavity wall. In dense shocks, dust grain 
temperatures are raised by collisional heating with hot 
gas particles and by absorption of photons as the com- 
pressed gas cools, resulting in thermal desorption and 
sputtering of interstellar ices. In particular, for a nondis- 
sociative shock with a velocity of v s = 10 km s" 1 , dust 
grains are expect ed to be heated to temperatu res on the 
order of 10 2 K (jNeufeld fc Hollenbachl Il99l . which is 
more than sufficient to allow for a condensation sequence 
to occur behind the shock. Therefore, similar to the ac- 
cretion outburst scenario, we propose that nearly pure 
CO2 ice mantles are formed in post-shocked gas, again in 
a condensation sequence, which the line-of-sight must be 
dominated by shock-processed material near the cavity 
wall. For other orientations, the line-of-sight will be dom- 
inated by irregularly-shaped icy grain mantles. Toward 
HOPS-68, evidence for the presence of shock-processing 



may be inferred from the detection of [S I], [Fc II], and 
[Si II] forbidden emission at 25.25, 25.98, and 34.82 fxm, 
respectively (D. M. Watson et al. 2013, in preparation), 
as well as gaseous absorption from C2H2, HCN, and CO2 
at 13.71, 14.05, and 14.97 /im, respectively. 

We note that the formation of interstel lar ices in a 
post-s hocked gas was originally proposed bv lBergin et al.l 
i|1998l 119991) . However, in their scenario, CO2 ice is 
formed by virtue of gas-phase reactions (CO + OH — > 
CO2 + H) and subsequent freeze-out. The formation 
process is largely governed by the photodissociation rate 
of H2O molecules by cosmic-ray induced photons, and 
is expected to occur over a timescale of t « 10 5 yr at 
the assumed density of nn 2 = 10 5 cm" 3 ; however, the 
timescale may be shortened with higher gas densities or 
larger cosmic-ray ionization rates. Nonetheless, this long 
timescale for CO2 ice mantle formation would allow for 
the re-development of irregularly-shaped grains, through 
the aggregation of small regular particles. In contrast, 
the detection of spherical icy grains implies that ice man- 
tle formation would have occurred more recently, over 
much shorter timescales (t ~ O(10 2 ) yr) and at higher 
densities (tih 2 ~ 10 7 cm~ 3 at r — 100 AU), by evapora- 
tion and re-condensation processes alone. 

In summary, we propose that the observed 15.2 /xm 
CO2 ice profile may result from the combination of two 
circumstances. First, the nearly complete absence of un- 
processed ices along the line-of-sight is due to the highly 
flattened envelope of HOPS-68, which lacks cold absorb- 
ing material in its outer envelope, and possesses a large 
fraction of material within its inner (10 AU) envelope 
region. Second, an energetic event led to the evapora- 
tion of inner envelope ices, followed by cooling and re- 
condensation, explaining the sequestration of CO2 in a 
hydrogen-poor mixture and the spherical shape of the icy 
grains. Regardless of which energetic process is respon- 
sible for producing the spherical mantles, we note that 
both proposed mechanisms require fortuitous timing and 
orientations, and hence are consistent with the observed 
rarity of spherical icy grains toward protostars. 

The notion of ices forming in a cooling gas, following 
an energetic event, makes several predictions. Water ice 
condensed at high temperatures will arrange itself into 
a crystalline matrix, predicting that the H2O ice seen 
toward HOPS-68 should be predominantly crystalline. 
Evidence for the presence of c rystalline H2O ice was pre- 
dicted, bv lPoteet et~aTl (|2011l) . from the 12 fj.m libration 
mode. This prediction can be further tested by observing 
the profile of the 3.08 fim H2O stretching mode. Carbon 
monoxide ice, observed at 4.67 /im, should either be ab- 
sent if the temperature of the absorbing material is still 
above 20 K, or should show evidence of spherical grain 
shapes, in shar p contrast to that observ ed toward typ- 
ical protostars (|Pontoppidan et all [20031 . Finally, the 
12 C0 2 and 13 C0 2 stretching modes at 4.27 and 4.38 /im, 
respectively, should be consistent with spherical grains 
and a nearly pure CO2 mixture. Using the best-fit mod- 
els from Sections 13.31 and 13. A\ predictions for the 12 C02 
and 13 C02 stretching modes are presented in Figure |U 
While HOPS-68 is too faint at 3—5 /im for current spec- 
trometers, these tests will likely be within reach of the 
James Webb Space Telescope. 
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Figure 5. The background- and baseline-corrected absorption spectra of pure CO2 ice. The (i> 3 ) 12 C0 2 and (S3) 13 C0 2 stretching (2343 
and 2283 cm -1 , respectively) and (P2) 12 C02 bending (660 and 655 cm -1 ) modes are shown at laboratory temperatures of T = 15 (black 
line), 30 (red line), 45 (blue line), 60 (green line), and 75 K (magenta line). Gas- phase CO2 at 668 cm -1 is a background-subtraction 
artifact, arising from cither under- or over-subtraction of background CO2 gas. 
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APPENDIX 

PURE C0 2 ICE LABORATORY SPECTRA 

The infrared absorption bands o f CO2 have been exten - 
sively studied in th e l abora tory bv lHudeins et al.l (119931). 
lEhrenfreund et all (|1997l ), and iBaratta fc Palumbol 
( 1998f ). Collectively, these studies present absorption 
spectra and optical constants for pure CO2 at labora- 
tory temperatures of T = 10, 12, 30, 50, 70, and 80 
K. These experiments are conducted under high- vacuum 
conditions at a spectral resolution of 1—2 cm -1 , too low 
to fully resolve the narrow structure of the Davydov split. 
More recently, CO2 experiments with a spectral reso- 
lution of 0.5 cm -1 have been c onducted at laboratory 
temp eratures of T = 15—90 K (jvan Broekhuizen et al.l 
120061 ). but at relatively low signal-to-noise. 

In this paper, we deploy new high-quality in- 
frared absorption spectra and optical constants from a 
temperature-series laboratory study of pure, crystalline 
CO2 ice. The experiments were conducted at an im- 
proved spectral resolution of 0.1 cm" 1 and at higher 



signal-to-noise with respect to previous laboratory stud- 
ies. 

Experimental Procedure 

Fo llowing the pro c edures of iGerakines et al.l (|1995[ ) 
and iBouwman et al.l (|2007f ). the experiments are per- 
formed in a high vacuum (HV) chamber with a base 
pressure of P « 5 x 10 -7 mbar at room temperature. The 
chamber houses a Csl (Caesium Iodide) substrate that is 
cooled down to 15 K by a closed cycle helium cryostat 
(ADP DE-202). The substrate temperature is controlled 
with a resistive heater element and a silicon diode sensor 
using an external control unit (LakeShore 330). 

A sample of CO2 (Praxair, 99.998 % purity) is intro- 
duced into the system from a gas bulb at 10.0 mbar pre- 
pared in a separate vacuum manifold (base pressure of P 
sw 10~ 5 mbar). Pure CO2 ices are grown onto the sub- 
strate at 15 K via effusive dosing along the surface nor- 
mal at a deposition rate of ^10 15 molecules cm" 2 s" 1 . 
The ices are heated to temperatures of T = 15, 30, 45, 
60, and 75 K at a rate of 2 K min -1 and allowed to re- 
lax at each temperature for five minutes before recording 
the absorption spectra. A Fourier Transform InfraRed 
(FTIR) spectrometer (Varian 670-IR) is used to record 
the ice spectra in transmission mode over the wavenum- 
ber range of v = 4000-400 cm" 1 (A = 2.5-25 fim) with 
a spectral resolution of 0.1 cm -1 , averaging a total of 256 
scans to increase the signal-to-noise ratio. Background 
spectra are acquired at 15 K prior to deposition and sub- 
tracted from the recorded ice spectra. 

Experimental Results 

To facilitate a comparison with astronomical spec- 
tra, a baseline correction is applied to all labora- 
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Figure 6. The &2 bending mode spectrum of pure CO2 ice (black 
line) compared to its corresponding "smoothed" spectrum (blue 
line) at T = 45 K. A least-squares fit is performed using a super- 
position of Gaussians (green lines) to eliminate noise and gas-phase 
CO2 absorption at 668 cm -1 . The resultant fit residuals (magenta 
line) are shown for comparison. 



tory absorption spectra by fitting a third-order poly- 
nomial over the wavenumber regions of 4000—3800, 
2900-2700, 2270-2250, 950-700, and 645-640 cm" 1 . 
The background- and baseline-corrected absorption spec- 
tra for pure CO2 ice are shown in Figure [5] for the (£3) 
12 CO and (£3) 13 CO stretching and (z> 2 ) C0 2 bending 
vibrational modes. 

The spectra do not suffer from any significant contam- 
ination by background gases of H 2 within the vacuum 
system. However, weak traces of gaseous C0 2 , near 668 
cm -1 (A = 14.97 /xm), are present in some absorption 
spectra, and are the result from either an under- or over- 
subtraction of background C0 2 gas. To reduce the noise 
and eliminate the gas-phase C0 2 absorption, a nonlinear 
least-squares minimization is simultaneously performed 
to the £3 12 CO, £3 13 CO, and £ 2 C0 2 vibrational modes 
using a superposition of Gaussians. Figure [5] illustrates 
the 'smoothing' technique for the £ 2 bending mode spec- 
trum at T = 45 K. The measured peak positions (D) and 
widths (FWHM; AD) of the z> 2 bending mode are listed 
in Table [3] The position of the absorption peaks near 
654.5 and 659.7 cm -1 are temperature-independent and 
remain constant to within 0.1 cm -1 , while their FWHMs 
decrease by a factor of two with increasing temperature 
from T — 15 to 75 K. These result s are in agreement with 
the 0. 5 cm -1 resolution study bv Ivan Broekhuizen et al~l 

(pool . 

Optical Constants 

The complex index of refraction (n + ik) is determined 
using a Kramers- Kronig anal ysis, following s i milar m eth- 
ods previously descri bed in iHudgins et al.l (|1993l ) and 
iMastrapa etaLl (|2009D . 

Briefly, the imaginary part of the complex refractive 
index, k, is calculated from the Beer-Lambert law: 



a(v) = knkv = - ln(10) A 9 , 



(Al) 



Table 3 

Band Positions and FWHMs of the V2 CO2 Bending Mode 



T 


V 


FWHM, Au 


Ice Thickness, d 


(K) 


(cm- 1 ) 


(cm" 1 ) 


(ML) a 


(jjeao) 


Laboratory Spectrum 


15 


654.5/659.7 


2.1/5.3 


366 


0.203 


30 


654.5/659.8 


2.0/4.3 


367 


0.203 


45 


654.5/659.7 


1.8/3.6 


410 


0.227 


60 


654.6/659.7 


1.5/2.9 


405 


0.224 


75 


654.7/659.7 


1.1/2.3 


355 


0.197 


CDE Model 


15 


654.9/660.6 


2.1/7.7 






30 


655.0/660.7 


2.1/6.6 






45 


655.0/660.7 


1.9/5.3 






60 


655.1/661.0 


1.7/4.6 






75 


655.3/661.0 


1.6/4.3 






Homogeneous Spheres Model 


15 


655.6/662.6 


1.4/3.6 






30 


655.6/662.5 


1.3/2.4 






45 


655.6/662.4 


1.1/0.9 






60 


655.6/662.3 


0.9/0.5 






75 


655.8/662.4 


0.6/0.3 






Ice-coated Spheres Model b 


15 


655.0/660.4/666.4 


2.0/13.7 






30 


655.0/660.4/665.8 


1.9/12.8 






45 


655.0/660.4/664.8 


1.7/ 9.6 






60 


655.1/660.5/664.4 


1.4/ 6.3 






75 


655.2/660.6/664.4 


0.7/ 4.8 







Note. — The band position and FWHM arc listed for the 
multiple peaks within the CO2 bending mode. 

a 1 mono-layer (ML) £j 10 15 molecules cm -2 . 
b For broad profiles exhibiting blended peaks, the FWHM 
value represents the overall width of the profile containing 
multiple peaks. 



where a is the Lambert absorption coefficient, d is the 
thickness of the ice, and is the absorbance of the ice 
at wavenumber v. The thickness of the ice samples were 
estimated using the known absorption band strength of 
pure C0 2 ice. In this technique, the column density, N 
(in molecules cm" 2 ), of the sample was determined using 



N 



A 



(A2) 



where tz = \n(10)Ap is the optical depth and A is the 
intrinsic band strength of the vi bending mode as deter- 
mined from previous laboratory studies (A = 1.1 xl0~ 17 
cm molecule -1 ; iGerakines et al.lll995[ ). The ice thick- 
ness is then inferred from the derived column density 
assuming a monolayer (ML) surface coverage of ^10 15 
molecules cm -2 . The final ice thickness is tabulated in 
Tableland ranges in value from 355—410 ML. Assum- 
ing a monolay er thickness of 5.54 A, th e lattice constant 
for solid C0 2 (|Keesom fc Ko hlcr 1934]), these values cor- 
respond to a physical thickness of ^0.20—0.23 /um. 

Given the absorption coefficient, the real part of the 
complex refractive index, n, is then calculated from the 
Kramcrs-Kronig dispersion relationship: 



n{y) = n + 3? 

Z7T 



4000 



100 



a(v>) 



dv' 



(A3) 
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Figure 7. The optical constants (n and k) for the V2 bending 
mode of pure CO2 ice at laboratory temperatures of T = 15 (black 
line), 30 (red line), 45 (blue line), 60 (green line), and 75 K (ma- 
genta line). 

where 2? denotes the Cauchy principal value of the 
integral and no is the refractive index for the high- 
wavenumber end of the infrared spectrum (y > 4000 
cm -1 ). Because the visible (A = 0.633 fim) refractive in- 
dex of pure CO2 ice is temperature-dependent, wc adopt 
values of no ~ 1.23, 1.26, 1.29, and 1.34 for the laboratory 
temperatures of T = 15, 30, 45, and 60 K, respectively 
(jSatorre et al.ll2008h . For the absorption sp ectrum at T 
= 75 K, we use a value of no ~ 1.44 from iSeiber et all 
(|1971h to ensure that n is positive for all wavenumbers. 
Optical constants for the v-i bending mode of pure CO2 
are presented in Figure [7J 

While our study is not intended to address the differ- 
ences between the CO2 optical constants derived in this 
work and those found among the literature, we do note, 
however, that any discrepancies are most likely due to 
differences in ice sample thickness and the no refractive 
index for the high-wavenumbcr region [v > 4000 cm -1 ). 
For a thorough discussion on the discrepancies between 
CO2 optical const ants derived from d i fferent studies, we 
refer the reader to lEhrenfreund et all ()1997l ) . 

Grain Size and Shape Corrections 

The profiles of strong infrared resonance features, in- 
cluding those of ices, depend significantly on th e shape 
and s i ze of the absor bing and scattering particles (jTielensI 
fl^9ltlDartoi^[27)06l ). For pure CO2 ice, grain size and 
shape corrections must be applied to the laboratory ab- 
sorption spectra of thin films to allow for an accurate 
comparison to astronomical spectra. While irregularly- 
shaped grain models have been largely successful at 
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Figure 8. The absorption cross section per unit volume for the 
i>2 bending mode of pure CO2 ice. Grain shape models in the 
Rayleigh limit are shown for a continuous distribution of ellipsoids 
(CDE; top panel), homogeneous spheres (middle panel), and ice- 
coated silicate spheres (bottom panel) at laboratory temperatures 
of T = 15 (black line), 30 (red line), 45 (blue line), 60 (green line), 
and 75 K (magenta line). 

simulating astronomical spectra of CO and C O2 ices 
(jPontoppidan et all 2003; Gcrak ines et al.llT999[ ). we also 
consider here the shape effects of spherical grain models. 

Using the derived optical constants of pure C Q2, stan- 
dard formulae from IBohren fc Huffman! (fl983h were ap- 
plied in the Rayleigh limit to calculate the absorp- 
tion cross section per unit volume, C a b s / V, for a con- 
tinuous distribution of ellipsoids (CDE; all shapes are 
equally probable), homogeneous spheres, and ice-coated 
amorphous silicate spheres. For the ice-coated spheres 
mod el, we adopt amorph ous silicate optical constants 
from IDraine k, Leel (|1984l ) and assume that half of the 
total particle volume is occupied by the core. Note, how- 
ever, for thick ice mantles or silicate cores that occupy < 
10% of the total particle volume, the ice-coated spheres 
profile approximates that of pure homogeneous spheres. 

The absorption cross sections for the bending mode 
of CO2 are presented in Figure [8] for the different grain 
shape models. The CDE model produces a broad double- 
peaked absorption profile that remains mostly constant 
with increasing temperature from T = 15 to 75 K; the 
peak positions exhibit a blueshift of 0.4 cm" 1 at T 
= 75 K. Conversely, a narrow double-peaked structure 
is generated by the homogenous spheres model, whose 
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Figure 9. The normalized optical depth profiles for the V2 bend- 
ing mode of pure CO2 ice. The uncorrected laboratory spectrum at 
T = 15 K (black line) is compared to three standard grain shape 
models in the Raylcigh limit: a continuous distribution of ellip- 
soids (CDE; blue line), homogeneous spheres (magenta line), and 
ice-coated amorphous silicate spheres (green line). 

FWHMs become narrower with increasing temperature; 
the FWHM of the 662.5 cm -1 peak decreases by a fac- 
tor of 12 from T = 15 to 75 K. In agreement with 
lEhrenfreund et all ()1997l ). we find that the ice-coated sil- 
icate spheres model gives rise to a broad triple-peaked 
profile. 

A comparison between the different grain shape mod- 
els and the laboratory spectrum of pure CO2 at T = 15 
K is presented in Figure |H1 For all models, the peak po- 
sitions and FWHMs deviate from those of the original 
laboratory spectrum. For the CDE model, we find that 
the profile becomes blucshiftcd and broadened by ~1 and 
~2.5 cm -1 , respectively. For the homogeneous spheres 
model, a large blueshift of ^3 cm -1 and a narrowing of 
~2 cm -1 are observed compared to the laboratory spec- 
trum. For a full comparison of the spectra, the reader is 
referred to Table [3j 
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